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Abstract 

We present here a new shear free model for the gravitational collapse of a spherically symmetric 
charged body. We propose a dissipative contraction with radiation emitted outwards represented by 
the Vaidya-Reissner-Nordstrom metric. The Einstein field equations, using the junction conditions 
and an ansatz, are integrated numerically. A check of the energy conditions is also performed. We 
obtain that the charge delays the Reissner-Nordstrom black hole formation and it can even prevent 
the collapse. 
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I. INTRODUCTION 

The equilibrium of the stars depends on the balance of two conflicting effects. The gravi- 
tational force attracting the material of the star toward the center, and the opposite internal 
thermal pressure provided by nuclear fusion of the elements in the stellar interior. When 
these reactions cease and no other source of pressure acts significantly, like the non-thermal 
degeneracy pressure, this balance is broken and a massive star undergoes the continuous 
catastrophic contraction, the gravitational collapse. The remnant product of this process 
can be a black hole. 

The General Theory of Relativity furnishes the foundation for analyzing the effects of the 
huge gravitational field involved in this process. In this context, the first idealized model 
for the gravitational collapse was proposed by Oppenheimer and Snyder in 1939. It 
consists on the dynamics of a simple pressureless fluid (dust), a set of no interacting parti- 
cles subjected only to the action of the gravitational field generated by their own masses. 
In consequence, the particles travel along timelike geodesies towards the center of the con- 
figuration, resulting in the formation of a Schwarzschild black hole. Thus, an interesting 
step forward is to verify the consequences in considering more relevant fluids, for instance, 
endowed with pressure, viscosity and/or electric charge. Since this landmark paper, many 
attempts aiming a more realistic description of this phenomenon have been taken place, but 
this usually impose a more complete set of equations, and the non linearity coming from 
the own structure of the field equations often compel us to the application of numerical 
procedures. 

An interesting issue lies on the fact that, if a star could hold a non null amount of net 
electric charge, the contraction of such an object could give rise to the Reissner-Nordstrom 
black hole. However, the existence of considerable fraction of net charge in the stellar 
..no. ha. been fo. Ion. .e.ected. A p.oce. proposed Homeland Q .... .Ka. .Ke 
free electrons of the strongly ionized gas that composes the star can be ejected due to its 
large thermal velocities. As the balance of forces is established, the process is interrupted 
and an object containing only 100 Coulomb per solar mass remains. Similar arguments 
supported by Glendenning |3[ reaffirm the neutrality of the stars. However, these objections 
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are grounded in the Newtonian context and, as it was already advertised j4|, do not apply 
when extremely huge gravitational fields are in play. Indeed, in the scenario of the Theory 
of the General Relativity, a relativistic star could hold greater portions of charge keeping 
itself stable (jj], jsl), albeit no mechanism producing charge asymmetry is known. A work 
in this direction can be found in jo]. Nevertheless, these issues seem far from be concluded, 
making it important to draw attention to the charged stellar models. 

There exist several relativistic static interior charged fiuid solutions connected to the 
exterior Reissner-Nordstrom spacetime^ A good classification and description of most of 
them can be found in Ivanov, 2002 [7]. Among these solutions, the one developed by 
Cooperstock and De La Cruz 8| will be useful for our investigation. It represents a charged 
dust fiuid, where the Tqq component of the stress-energy tensor in the comoving frame 
- composed by the matter and the interior electric field - is a constant, giving rise to a 
Schwarzschild-like interior solution. 

In spite of the several static charged interior metrics, dynamic solutions are rare. In this 
sense, an interesting work performed by Bekenstein in 1971 [9] has indicated the dual role 
played by the charge. Depending on the charge-mass ratio, the electric field can contribute 
to the support or, in opposition, to the speeding up of the structure's collapse (charge 
regeneration). The proof of the Birkhoff theorem and the hydrostatic equilibrium equation 
for the charged case were also achieved. 

The Vaidya-Reissner-Nordstrom metric, generalization of the Vaidya's metric 10| in the 
presence of electromagnetic field, enables the geometrical description of the radiation field 
surrounding the star, allowing the development of several models of radiative gravitational 
collapse of charged star s |lll -ll4| . Among these models, the non adiabatic shear free model of 
de Oliveira and Santos [ill] is very similar to the one proposed in the present paper, despite 
the fact that the field equations have not been integrated before and of the absence of pressure 
anisotropy. Some years earlier, the same authors and Kolassis [15[ studied the uncharged 
case by means of the separation of the metric functions of the general spherically symmetric 
isotropic interior metric into distinct functions of the radial and temporal coordinates, a 



proposition already known from a preceding investigation [16|. Years later these models 
have been generalized by Chan for the shearing model 17|]. This simplification allows the 
complete integration of the Einstein field equations, provided a radial function representing 



the static pre-collapse configuration is given. Here we use this same assumption, besides the 
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selection of the referred Cooperstock and De La Cruz solution js] for the the radial function 
that represents the initially static situation. We are going to see that this ansatz forces the 
physical variables out of equilibrium to have the same radial behavior as if they were in 



equilibrium. This is called in the literature as the post-quasi-static approximation |18N20|. 

In studying a static relativistic stellar model, Bowers & Liang I2I] pointed to the relevant 
effects played by the anisotropy of pressure over the maximum mass for a stable system 
and over the surface redshift. among other factors, pressure anisotropy is induced by the 
electric charge 22(|, by the shear movement 17|,|23| and, moreover, it is enlarged by the shear 



viscosity 



24l | (see 25| for a detailed review article about the occurrence of local pressure 



anisotropy in the study of self-gravitating systems). The shear movement, on the other hand. 



affects the luminosity and total collapse time of the contracting systems [26|,l27|, justifying its 
inclusion in several works. However, in an unsuccessful attempt to obtain an exact solutions 
of the Einstein's equations, in this paper we forsake the general case of shearing collapse, 
although, due to others sources, the pressure anisotropy is taken into account. 

Here we have to remind important consequences of the shear-free condition already stud- 
ied. It was shown by Herrera and Santos 28|, l29| that, in the quasi-static approximation 
and for the non-dissipative situation, the shear-free flow is equivalent to the homologous 
contraction in the Newtonian limit (details about homologous contraction, homology rela- 



30|). Moreover, considering 



tions and its application in the astrophysics can be found in 
dissipation in both regimes, diffusion and streaming out, the authors verified the need to 
impose homology conditions on the temperature and emission rate to keep the homologous 



evolution [28|. In addition, it was checked that the requirement of homogeneous expansion 
implies homology conditions. Another important aspect is that the shear-free condition is 
not assured all along the collapse of spheres whether the dissipative processes, the density 
inhomogeneity and the pressure anisotropy are taking in account 31|- Hence, this follows, 
leastwise in the study of geodesic motion of the fluid. 

Current literature is comprised by some comprehensive models of charged collapsing 
radiating spheres with shear motion. A good example is the work developed by Medina and 



collaborators 



12| with the aid of a method developed by Herrera et. al. [32]. Following the 



procedure, a system of first order differential equations evaluated at the surface is obtained 
and integrated. The .a^e n,ethod was appUed by Resale, et. al. R in .ecent pape. 

Treating the same subject, but using the approach of Misner and Sharp [33], Di Frisco and 
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collaborators 



34j analyzed the dynamics of the charged spheres. 



The stellar contraction must be a highly dissipative process due mainly to the interaction 
of the radiation and matter in the interior, such that, in the same way we are considering 
emission of energy across the stellar surface, an outward flux of energy through the body oc- 
curs. The energy transport is often considered in both regimes, diffusion and free-streaming. 
Once the studies of the radiation of the type-II supernova 1987A pointed to the domain of 
the former [sS^, this is the only one taken into account in this work. Details on how the 
energy flux happens and how it relates to the temperature are not the focus here. This 
could be performed, at first glance, with the usage of the relativistic thermodynamic theory 



developed by Eckart 



36(1 . But, although simple, this theory has an undesirable feature. 



the equations for the heat transfer and viscous tensions indicate instantaneous propagation 
of the perturbations, generating causality problems. This difficulty was overcome with the 
Israel- Stweart theory 37|], where the relaxation times for the dissipative quantities were in- 
cluded, resulting in hyperbolic equations for the perturbations (a complete discussion can be 
found in jsSi). The energy transport equations of this theory supplied the thermodynamic 
analysis in several works. In an interesting paper, Di Prisco and collaborators 39[ compared 
the luminosity profile of the collapsing spheres for the vanishing (Eckart) and non- vanishing 
(Israel-Stweart) relaxation times. The same idea was used in a different model by the same 
authors and some more collaborators [40] one year later. In this way, some effects of the 
pre-relaxation process were confirmed to be model-independent (more on the infiuence of the 
relaxation time in the outcome of the collapsing systems in 4l|)- In applying the transport 
equation to collapsing sphere models, the authors usually simplify the problem making null 
the constants coupling heat and viscosity. Few years ago, Herrera and collaborators j42| . 
using the the method by Misner & Sharp [33], studied the dissipative gravitational collapse 
when this simplification is set aside. In this way, the authors have appreciated the infiuence 
of these coupling constants over the effective inertial mass of the system. 

The aim of the present investigation is to find a dynamical solution for a charged sphere 
representing the dissipative collapse of the charged body, as well as to verify the role played 
by the electric charge in this process. The formation or avoidance of the black hole is 
discussed. The paper is organized as follows, in Section [III we set the metrics and the 
stress-energy tensor used, we show the field equations and obtain a set of equations by the 
application of the junction conditions. The solution of the system of equations is performed 
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in Section IIIII The equations for the static initial configuration are presented in Section IIVI 
The resuhs can be found in Section |V] and the agreement with the energy conditions are 
checked the Section |Vll Finally, we discuss the results in Section IVIIi 



II. EQUATIONS 

The system is composed by an interior spacetime (referred by a minus signal "-"), a 
comoving timelike hypersurface (referred by S) and an exterior spacetime representing a 
radiation field surrounding the body (referred by a plus sign "+"). The unknown inte- 
rior spacetime is given by the spatially isotropic spherically symmetric shear-free metric in 
comoving coordinates, 



ds^_ = -A\r, t)de + 52 (r, t) [dr'' + r''{d9'' + sin^ Odcj)'')] 



(1) 



The exterior spacetime is described by Vaidya-Reissner-Nordstrom metric, which repre- 
sents an outgoing radial flux of radiation around a charged spherically symmetric source of 
gravitational field, given by 



dsl 



1 - 



2m{v) ^ Q 



dv^ - 2dvdr + r^{de^ + sin^ 9d(j)^), 



(2) 



r 

where m{v) represents the mass, function of the retarded time v, and Q is the total amount 
of charge of the system inside the boundary surface S. 

The interior stress-energy tensor of a charged dissipative anisotropic fiuid is given by 



Z 



al3 



(/i + p±)UaUi3 + P±gal3 + (P " + 



An 



A" 



7(5 



(3) 



These quantities are the energy density /i, the radial pressure p, the tangential pressure 
p_L, the radial heat fiux g", -u" is the four- velocity and X" is a radial four- vector. The 
following relations hold: Uag" = 0,XqX" = = 0. The tensor F^j^ represents the 

electromagnetic field tensor. At last, the coupling constant in geometrized units is k = Svr 
(i.e., c = G = l). 

As we use comoving coordinates, we have 



(4) 
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and as the heat flux is radial 

= qSl (5) 

The Maxwell's equation are written as 

Fafi = 0/3,a - 0a,/3, (6) 

and 

F"^^ = Anr, (7) 

where (pa is the four-potential and J" is the four-current. It is assumed that the charge is at 
rest with respect to the coordinates of the metric ([1]), so we have no magnetic fleld present 
and, thus, we can write 

0. = (8) 

and 

r = au"". (9) 

Where $ is the scalar potential and a is the charge density. 
From equations ([8]) and ([6]) we obtain 

^01 = -i^io = (10) 
Substituting equations @ and f|TO|) into Maxwell's equations ([6]) and ([7]), we have 

$" + + 1^ + - V' = 47rAi?V, (11) 



A B r 



and 



where the dot represents differentiation with respect to t and the prime with respect to r. 
Integrating these equations, we obtain 

where l{r) is the radial distribution of charge given by 

r 

l{r) = 4vry" aB^r^dr. (14) 
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The integration from the center (r = 0) till the surface (r = r^) provides us the total amount 
of charge Q inside the body. 

The non-vanishing components of the field equations applied to the interior spacetime 



are very similar to the paper ll|], except by the presence - justified in the first section - of 
the anisotropy of pressure. These are obtained with the aid of the equations ([T]), OH]), (H]), 
(E]), resulting in 

2 
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B \B 



AB' 



3 1 -) =^,(A'fi + ——J-\, (15) 
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AB 
AB 



(16) 
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52 
A^ 



B" 


(By iB'_ 
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lA' 


B 


\B r B 


' A- 
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r A 




2^^ 

AB 



kB' 



871 B^r' 



(17) 



01 



B' B' B A' B , . 



(18) 



We can note that the equations presented in this section are a particular case of those 
present in the section II of the paper 3^. Here we are assuming no shearing motion and we 
neglect the shear viscosity as well as the dissipation in the free-streaming regime. 

We consider a spherical surface with its motion described by a timelike three-surface 
S, which splits spacetimes into interior and exterior manifolds. The connection between 
these geometries has to be continuous and smooth, assured by the junction conditions. 
Here we follow the approach given by 43|, 4^ and demand the continuity of the metric 



((is^)E = and of the the extrinsic curvature K^, = across the surface S 



(details on the junction conditions used here can be seen in the paper 11|). 
The condition of continuity of the metric leads us to 



dt 



A{rj:,t)-\ 
8 



(19) 



and 



2m dr 

1 + 2:r + ^ 



(20) 
(21) 



where r is a proper time defined on the comoving hypersurface S. 

The non- vanishing extrinsic curvature components are K^^, Kf^ e K^^ 
Equating K^q and K'^q we have 



2m 



dr 
r + — r 

dr 



[r {rB' + B)], 



With the help of equations ( !T9|) . (!20|) . ( |2T|) . we can write equation (122!) as 



m 



2^2 



3d/2 



2fi 



+ 



2r5 



which is the total energy entrapped inside the surface S 
Equating and i^^, using equation (ITU]) , we have 



(22) 



(23) 



45|. 



d^v f dv\ ^ f dv\ 



m 



(ir^ \dT J \dT J r"^ c/r 



A' 



(24) 



Substituting equations ( 1T9|) . ( l20l) and ( l23i) into ( l22i) we can write 

-1 



dv 
dr 



{rB)' ^ {rBy 



(25) 



B A 

This is the gravitational redshift for the observer at rest at the infinity. Its divergence 
indicates the formation of an event horizon. This happens when the factor in parentheses 
goes to zero. 

Substituting the derivative of this last equation and equation fl23|) into identifying 
the geometrical terms with the field equations ( 1TB]) and (ITS]) , we have 



(26) 



This result is the same of the one obtained by de Oliveira et. al. [l5| and years later in 



the study of charged spheres 



Ul and (l3|. 



We can use the equations ( 123]) . ([19]) . ( 125]) and (|T6]) to obtain the total luminosity for an 
observer at rest at infinity 

2" 



dm 
dv 



K 

2 



pr^B"^ 



rB' rB 
1 + ^ + — r 



B 



A 



(27) 
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III. SOLUTION OF THE FIELD EQUATIONS 

In order to integrate the field equations, we resort to the same ansatz referred in the 
Section [H consisting in separate the metric functions A{r, t) and B{r, t) into functions of the 
coordinates r and t in the form 



where Ao{r) and B^ir) are solutions of a static charged dust fluid. 

We have chosen this separation of variables in the metric functions in such way that 
when f{t —7- — oo) — )■ 1, the metric functions represent the initial static solution. The 
decreasing (or increasing) of /(t) provides the decreasing (or increasing) of the physical 
radius (i? = rB{r,t)) of the body, describing the collapse (or expansion). The process ends 
when f{t ^ 0) ^ 0. 

Now the equations (fT5|) - (fT8|) can be written as 



A{r,t) = Ao{r), 



(28) 



and 



B{r,t)=Bo{r)f{t) 



(29) 




(30) 




(31) 




(32) 



Kq = 



2^0 / 
Ao'Bo'P' 



(33) 



where 




(34) 




(35) 
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and 



Bl 
B. 



Bo 



rBo Ao r Ao 



Substituting equations (13T]) . (!33|) into (l26l) and assuming also that poijY}) 
a second order differential equation for /(t), 

2/7 + /V'-2a/V-&(l-/') =0, 



where 



and 



5o 



(36) 

0, we obtain 
(37) 

(38) 

(39) 

except by 



This equation is similar to the one obtained by de Oliveira et al. in 1985 
the last term. The attempts to solve it analytically proved to be unfortunate, at last, the 
numerical solution was successful. In order to verify the electric field action, different values 
of total charge were considered. 



IV. MODEL OF THE INITIAL CONFIGURATION 



We consider that the system at the beginning of the collapse has a static configuration 
of a charged dust fluid described in Cooperstock e De La Cruz js]. The Schwarzschild-like 
interior solution is obtained considering 



rpO 
^0 



A + 



(40) 



l-K Svr ' 

where Tq is the mass-energy density component of the stress-energy tensor, A is the mass 
density, E is the electric field and a is a constant. 

The metric for this case, presented here in isotropic coordinates, is given by 



dsl = Al{f)de + Bl{f) [rff2 + f2(ci^2 + sin^ 



where 



2\2 



(41) 



(42) 
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and the constants (3 and a are given by 



a2C2(l +f2)2' 

f = r/C, 



(43) 
(44) 



,3 = 1- a^R^ + (1 _ a'RY', 



2 _ 2mo _ 

The quantities mo, -R = (t'e-Bo('^s)) and Q are the total mass, radius and charge of the 
object. The constant C comes from the integration of the original Schwarzschild coordi- 
nates transformed into the isotropic coordinates. It can be easily determined by junction 
conditions. 

The coefficients (l38l) and fl39l) become 

"'^^ (45) 



(1 + Vl-a2i?2)- 



. 2 ■ 



Substituting the functions Ao{f), Bq^t), their derivatives and a rewritten version of the 
equations (l37|l into the equations (130|) - ( 135|) . we obtain 
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/V /32/2(l + f2) 



r/ 



2^ ' 



2^4 



/2a 



up = — 



Kp± 



2^4 



/2a 



+ f2 



+ f2)' 



6 1- 



b 1 



-2a/ 



2a/ 



where 



2a2Cf / 



a2 (3 - f2) 



57r 



(46) 



(47) 



(48) 



(49) 



(50) 



Po{r) =P±o{r) = 0. 
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The two last equations correspond to the charged dust solution of the paper 
isotropic coordinates. For the charge distribution we have 

16f6 



but in 



(51) 



V. RESULTS 

The integration of the differential equation fl371) depends on the coefficients a and b 
(equation (l45l) ). related to the values of mass, radius and charge of the object before the 
beginning of the collapse. We have taken mo = 5Mq, R = 15 km and different values for 
the ratio charge-mass Q/rriQ ^ 0.00, Q/mo = 0.20,0.40,0.64. The result is shown in Figure 
[H from which we notice that the collapse is slower with the increasing of the charge. This 



result is in contrast to that obtained by Medina et al. 
authors. 



12|, pointed as unexpected by the 
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FIG. 1: Temporal behavior of the function f{t) for the models with Q/niQ ~ 0.0,0.20,0.40. The 
time is in seconds. 



In Section [TTl we declared that the formation of an event horizon occurs when the surface 
gravitational redshift goes to infinity, i.e., the term in the parentheses in the equation fl25|) 
goes to zero. Using equations ( 125]) . and the metric functions and into the 
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parentheses of that equation, we have 



1 + 



rB, 



_^ rBof 



1 



Bo Ao 
So, the event horizon arises when 

fbh = 



1 + 



R{1 + VI - a^R^) „• ' 
+ ^ hk 



0. 



(52) 



, 2 ■ 



(53) 



R{i + vi-am^y 

Note from Figure [2] that, when the amount of charge is not very large, the event horizon 
is always formed, however, for larger values, this condition is no longer satisfied. 



Q/mn = 0.630 



df/dt for Ihe AH formation = -2227 



Q/rrin = 0.631 



df/dt for ttie AH formation = -2233 



Q/m„ = 0.632 



df/dt for tfie AH formation = -2240 
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FIG. 2: Gravitational redshift condition for a star with initial mass 5 M© and initial radius R = 15 
km. For Q/rriQ > 0.6310, the function / no longer reaches the value for the event horizon formation. 



In the last stages of the contraction of the overly charged stars {Q/rriQ > 0.6310), the 
decreasing of the function f{t) decelerates and it reaches a minimum value (see the Figure 
|3]). Beyond this point, the function increases sharply and a bounce of the system seems to 
occur. This idea is false, since the luminosity becomes negative for positive / (see equation 
f l55|) below). This is physically inconsistent, and we have to stop the integration at this 
instant. The solution is just valid until this minimum point, where the luminosity, heat 
flux and the rate of decrease of the stellar radius (/) go to zero. Nevertheless, the mass 
is reduced along the contraction by the radiative emission, but some portion still lasts at 
the endpoint. The same is truth for the stellar radius. This suggests that an equilibrium 
situation could be the final fate of the evolution. 

The equations ( 123|) and (!27|) become 



m 



2/3'^ 



2Rf 



(54) 
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FIG. 3: Time evolution of f{t) for a highly charged star. 



L = -aR^f 



(55) 



Note that, before the collapse, when f{t) — t- 1 and / — ?■ 0, the equation fl54l) becomes 



mo 



2R' 



(56) 



as we expected from js]. 

The Figure m shows the temporal evolution of the mass and luminosity. We can note that 
models with larger values of charge favor higher emission of energy. However, in the case 
where we have no black hole formation, this situation is reversed, as we can see from the 
right bottom panel. 

The integrated physical quantities of the fluid are plotted in the figures O [6l [7] and HI 
Due to the fact that the charge distribution is an increasing function of r, for the charged 
case, the radial pressure in the outermost radial coordinates are greater than the pressure 
in the innermost. 



VI. ENERGY CONDITIONS FOR THE CHARGED FLUID 

In principle, any geometry might represent, via Einstein's equations, a valid solution of 
the theory of General Relativity, once nothing is inferred about the matter content of the 
source of the gravitational field, or, precisely, no restrictions are set on the components of 
the stress-energy tensor. On the other hand, the known forms of matter are endowed with 
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FIG. 4: Mass and luminosity for different charge/initial mass ratio models in geometric units. 
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FIG. 5: Temporal evolution of the density for an almost uncharged and a highly charged star in 
geometric units. 



certain physical features, like positivity of the energy density or like its domain over the 
principal stresses, for example. Desired physically meaningful solutions, these features are 
combined in a set of constraints, the energy conditions, guiding the acceptance or disposal 
of a particular solution of the field equations. 

In order to trace the energy conditions along the moment of contraction of the body, 
we follow the same procedure used in Kolassis, Santos & Tsoubelis ||46j. Here we skip the 
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FIG. 6: Temporal evolution of the pressure for an almost uncharged and a highly charged star in 
geometric units. 
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FIG. 7: Temporal evolution of the tangential pressure for an almost uncharged and a highly charged 
star in geometric units. 
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FIG. 8: Temporal evolution of the heat flux for an almost uncharged and a highly charged star in 
geometric units. 



mathematical step by step, warning that the similar approach is found in [231 . 14 7j . These 
conditions for the charged fluid in question are fulfilled if the following inequalities are 
satisfied: 

Mathematical condition 
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(i) Va>0 
Common to all the energy conditions 



1 P(r) f 1 



> 



Weak conditions: 



{in) 



Dominant conditions: 



1 '"('•) 



> 



{iv) 



1 P(r) f 1 



> 



Strong conditions: 



Where 



[VI] 



1 l\r) 



> 



1 P(r) 
A + 2p+ ^ > 0. 



A 



1 l^ir] 



8n r*B^_ 

The Figure M shows that all these inequalities hold all along the dynamics of the charged 
fluid. 
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FIG. 9: Energy conditions for the highly charged model. 
VII. DISCUSSION 

In this work wc have studied the contraction of a spherically symmetric charged body. 
In order to make it possible the integration of the field equations, we have to adopt some 
hypothesis about the interior spacetime. Firstly, its associated metric can be split into 
distinct functions of the radial and the temporal coordinates. Even though we can not 
guarantee the validity of this mathematical trick, it is a reasonable attempt to get round 
the underdetermination problem of the system of equations. This drives us to the second 
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assumption, the choice of the radial interior function. Among the uncountable possibilities 
of interior charged solutions, the Cooperstock and De La Cruz [s] solution was the selected 
one. Thus, the radial behavior of the physical variables of the system is an inheritance of 
such a particular choice. 

In view of the previous discussion, we have concluded that appreciable physical effects to 
the collapse process are important for values of the net electric charge above Q/itiq ~ 0.10 
{Q ~ 8.570 X 10^^ Coulomb). Furthermore, the Reissner- Nordstrom black hole is formed 
unless the total amount of charge is incredibly huge {Q ~ 5.408 x 10^° Coulomb for Q/rrtQ = 
0.631). 

The electric field delays the event horizon formation, in agreement with the non radiative 

n 

models of [5|, and can even prevent the complete contraction of the body. In this case, the 
solution seems to point to an equilibrium situation. It is even worthy to say that, for a 
model with ratio Q/tuq ^ 0.70, neither a minimal contraction takes place. 

We also have seen that, for the models that admit black hole formation {Q/rrio ^ 0.631), 
those richer in charge have greater transmition of energy to the exterior region before the 
appearance of the event horizon, making the remnant object less massive. The luminosity 
increases sharply and, all of a sudden, the star turns off. Otherwise, if it is charged enough 
in such a way to prevent the black hole formation, the time evolution of the luminosity has 
a maximum peak and a subsequent reduction to zero. In contrast, the event becomes less 
luminous with the increasing of the fraction Q/rrio- 

Although the model presented here is, for sure, too much idealized, it certainly represents 
an interesting dynamical solution of Einstein's equation, and it gives us important clues on 
the possibility of the Reissner-Nordstrom black hole formation. 
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